function matmemtst()
% matmemtst:  test the matrix memory operations...
%
dt=0.1;
t=0:dt:1000;
f = 3.25*10^2;
a=3*sin(2*pi*f*t);
b=gradient(a);
ao=a;
bo=b;
mem(0);
[c, d] = change(a, b);
mem(0);
ai = any((a-ao) ~= 0);
bi = any((b-d) ~= 0);
if(ai > 0)
   disp('a changed')
else
   disp('a unchanged')
end;
if(bi > 0)
   disp('b changed')
else
   disp('b unchanged')
end;



function [ao, bo] = change(a, b)
ao = [];
bo = [];
%ao(1)=-1000;
bo=b*1000;
mem(0);
return;
